x1=13;
y1=interp1(x,y,x1,'spline')

y1 =

   27.8725

xi=0:1/3600:24;
yi=interp1(x,y,xi,'spline');
plot(x,y,'*',xi,yi);
